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The recent result on the third order correction to the Higgs boson production through gluon 
fusion by Anastasiou et al. T| not only provides a precise prediction with reduced scale uncertain¬ 
ties for studying the Higgs boson properties but also establishes the reliability of the perturbative 
QCD. In this letter, we propose a novel approach to further reduce the uncertainty arising from the 
renormalization scale by systematically resumming the renormalization group (RG) accessible log¬ 
arithms to all orders in the strong coupling constant. Our numerical study based on this approach, 
demonstrates a significant improvement over the fixed order predictions. 
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The remarkable discovery of the Higgs boson with a 
mass of about 125 GeV by the ATLAS and CMS col¬ 
laborations [2] at the LHC has provided an important 
clue to understand the mechanism of spontaneous sym¬ 
metry breaking within the framework of the Standard 
Model (SM) of particle physics. The technological ad¬ 
vancements in experimental sectors augmented with the 
precise theoretical predictions, played crucial role in this 
distinctive discovery. But, with the new data to be avail¬ 
able soon at the upgraded LHC, minimizing the theoret¬ 
ical uncertainties will be of paramount importance. The 
pursuit of the precision studies in the Higgs boson pro¬ 
duction has been a consistent pioneer in advancing the 
perturbative QCD. It is worth recognizing the fact that 
the fixed order [3] as well as threshold resummed [I] pre¬ 
dictions in perturbative QCD along with the electroweak 
effects [5] played an important role not only in the ex¬ 
clusion of wide range of the Higgs boson masses but also 
to establish that the discovered boson is almost consis¬ 
tent with that of the SM. Recent computation [5] of the 
complete threshold corrections at next-to-next-to-next- 
to leading order (N 3 LO) including the (5(1 — z ) part has 
marked a milestone. Owing to the universality of the soft 
emissions, this result was followed by various new results 
Hi for QCD processes at N 3 LO in the threshold approx¬ 
imation. Very recently a state-of-the-art computation 
J] has been performed by Anastasiou et al. to accom¬ 
plish the complete N 3 LO perturbative QCD correction to 
the inclusive Higgs boson production in the gluon fusion 
channel. This N 3 LO corrected result not only demon¬ 
strates the reliability of the perturbation theory through 
the moderate correction, but also reduces uncertainties 
significantly resulting from renormalization (fin) and fac- 
torization scales in the range /i £ where 

toh is the mass of the Higgs boson. Up to next-to-next-to 
leading order (NNLO), it was demonstrated in [8] that 
there was a significant increase in scale uncertainties if 
we increase the range. We also observe a similar pattern 
even at N 3 LO level for the fj,R variation. This is because 
of the presence of large logarithms of the scale at every 


order. Resumming such logarithms could often improve 
the scenario. In this letter, we use RG invariance of the 
Higgs boson production cross section to systematically 
resum these large logarithms to all orders in perturba¬ 
tion theory and show substantial reduction in the scale 
uncertainties over the fixed order predictions. In [9], for 
the Higgs boson production, it was shown that the large 
corrections of the form (Cattcl s ) 2 resulting from analyt¬ 
ical continuation of the form factors to time like regions 
can be successfully resummed to all orders using RG, 
giving rise to reliable predictions for K factor. Using 
effective held theory approach, the authors of [10] have 
shown the role of RG in improving the theoretical pre¬ 
dictions. Our approach, while uses same RG invariance, 
differs from theirs in treating the expansion parameter 
in a systematic manner as it will be demonstrated in the 
following. 

The inclusive hadronic cross section (<r H (s,m^)) for 
the Higgs boson production is related to the partonic 
cross-section as 
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where, f a {xi,p 2 F ) and fb(x 2 ,p 2 F ) are the parton distri¬ 
bution functions (PDFs), renormalized at /ip-, of the ini¬ 
tial state partons a and b with momentum fractions x± 
and x 2 , respectively and r = m^/s with y/s being the 
hadronic center of mass energy. a s = a s /4-7r with a s 
being strong coupling constant, cr° is an overall factor 
describing the effective interaction between gluons and 
the Higgs boson at lowest order and Cr is the Wilson co¬ 
efficient. Expressing a K (s,m'^ i ) = a 2 (^) a(s, /rfj), 
and using the RG invariance of (t h (s,to^), namely 

V%dk aH = °> we find 
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where, /3(a s (/r 2 )) = ^ 2 ^a s {^ 2 ) = - E^o A «i +2 (^ 2 )- 
Considering /uo as the central scale and using naive evo¬ 
lution of a s , Eq. [2] can be solved order by order to obtain 
the following perturbative expansion of 
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where = In (^§Y The coefficients of logarithms at 

each order in a s , lZ n ,k(0 < k < n) are governed by the 
RG evolution and can be expressed in terms of the lower 
order ones, lZ n -i , 0 , through 
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The coefficient of the highest logarithms at n th order in 
a s grows as (n + l)a"/?o 72o,o which often can give rise 
to potentially large contributions and can make the fixed 
order predictions unreliable. The RG invariance can be 
used to resum such contributions to all orders. To achieve 
this task, we extend the approach of m to the case of 


scattering cross sections in hadron collisions. We rewrite 
Eq. [3] as 
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so that <rE resums a s (^ R )L R to all orders. The closed 
form of can be obtained using RG invariance. The 
recursion relations (Eq. [4]) which follow from the RG in¬ 
variance, can be used to show that E"”* satisfies the fol¬ 
lowing first-order differential equations 
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where 0 m _i is Heaviside Theta function, 77 ^ = /3i/ /3 q and 
uj = 1 — /3oa s (fj, 2 R )L R . Upon solving the above equations 
recursively, we obtain a m for all m. In Eq. [7] we present 
them up to m = 4. 
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Alternatively <7^^ can be computed from Eq. fusing RG improved solution for a Sl given in Eq. [sj which implicitly 
resums the large logarithmic contributions to all orders in the perturbation theory. 
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In the above Eq.(|8|, the terms up to a 4 is already known, 
[12] and a 3 term is obtained for the first time. In the 
following, we study the numerical impact of fixed or¬ 
der (FO) as well as RG improved resummed (RESUM) 
cross sections up to N 3 LO in QCD for the Higgs bo¬ 
son production through gluon fusion at the LHC. We 
have used an in-house Fortran code to do this. We 
set 7 t 0 = hf = ttih = 125 GeV throughout and 
use MSTW2008nnlo [T3] parton distribution functions 
with the corresponding strong coupling constant from 
LHAPDF [H], a s (mz ) = 0.11707. At LO, the exact 
top and bottom quark mass effects are included through 
cr° in Eq. 0 Finite quark mass effects at NLO are taken 
into account using iHixs at hr = hf- At NNLO and 
N 3 LO, we use effective theory predictions in the large 
top quark mass limit. We first obtain Lr independent 
terms namely -Ro.Oi-Ri.d an< i 7?2,o by setting hr = ulr in 
our code whereas R 3,0 is extracted from the recent result 
for N 3 LO cross section given in (1. for the same choice 
of hr = hf = ttih- These Ri^o, (i = 0,1,2, 3) thus ob¬ 
tained at hf = 777h with MSTW2008nnlo are the only 
required ingredients to study the hr dependence of both 
the FO (Eq. [3]) and the RESUM (Eq. 0 cross sections 
up to N 3 LO in QCD. Note that the coefficients of all the 
Lr s in Eq. B can be obtained using the recursion rela¬ 
tions (Eq. [4]) . As it was demonstrated in [Tj, inclusion of 
N 3 LO corrections makes the hr sensitivity of the cross 
section milder compared to NNLO corrected results when 
the hr is taken to be closer to mn, say between tuh/4 
and 2t77h- On other hand, if we decrease hr. below 777.h/4, 
the contributions from Lr increase substantially surpass¬ 
ing the scale independent ones giving rise to potentially 
large scale uncertainties. This happens at every order 
in perturbation theory and the renormalization scale at 
which this happens, increases with the order. In Fig. |T} 
we quantify this up to N 3 LO for FO. 

In the FO results (Eq. [ 3 ]), the dependence on hr enters 
through the evolution of a s (/x 2 j) as well as the perturba¬ 
tive corrections that are polynomials in Lr of the order 
k < n consistent with RG invariance. As hr decreases, 
the coupling constant as well as the magnitude of Lr 





FIG. 1: Hr dependence of the LO and higher order correc¬ 
tions (FO) for LHC13, keeping hf = rn h fixed. 


will increase, consequently, for hr much less than m h, 
the contributions of the kind )L k R can become large 
enough to make the hr dependent terms even negative. 
Moreover, at higher orders, the contribution of polyno¬ 
mial in Lr need not be monotonic, instead it can change 
its sign with decreasing hr as seen in fig(T] With in¬ 
crease in Lr, the presence of the terms (/3oa s Lfl) fc makes 
the truncation of the perturbation series unreliable. The 
solution, proposed in this letter, resulting from RG im¬ 
proved resummation of those terms that spoil the pertur¬ 
bation series, shows an impressive improvement at every 
order. In figj2j we show both the FO and the RESUM 
cross sections up to N 3 LO for LHC13 by varying hr in 
the range [O.Dtih, IOtor] and keeping hf = tuh fixed. 
For hr < as discussed before, the large contribu¬ 
tions from Lr make the FO QCD corrections flip the 
sign and hence the cross sections take a downturn be¬ 
low certain hr- This phenomenon can foremost be seen 
for cross section at higher orders, e.g., for hf = 7tih the 
N 3 LO cross section starts declining below hr = 0.5mH, 
followed by NNLO cross section at hr = 0.2t71h and so 
on. For hf — 2toh also a similar pattern can be seen. 
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For larger values of /ir > mu, however, a a (/x|j) falls down 
suppressing the logarithmic contributions and hence the 
cross sections will decrease monotonically. We have also 
plotted the RESUM cross sections at various orders in 
Fig- i as a function of hr- We find that the predictions 
from the RESUM cross sections are more stable com¬ 
pared to the FO ones over a wide range of hr demon¬ 
strating the power and the reliability of resummation. 



FIG. 2: hr dependence of both the fixed order and resummed 
cross sections up to N 3 LO. 



LO 

NLO 

NNLO 

N 3 LO 

FO (%) 

167.26 

143.40 

54.99 

27.01 

RESUM (%) 

6.11 

5.47 

3.39 

1.23 


TABLE I: Percentage of maximum uncertainty for hr varia¬ 
tion in the range [O.lmn, lOmn] up to N 3 LO (see text). 

In Table [TJ we show the maximum percentage of uncer¬ 
tainty in the cross sections up to N 3 LO for hr variation 
in the range [O.lmn, lOmn]- Here, at N 3 LO, the hr un¬ 
certainty is maximum for hr between about O.lmn and 
0.5?7in whereas at NNLO, the maximum uncertainty is 
for /lr between about 0.2mn and lOmn- We notice that 
the scale uncertainties in both FO and RESUM cross sec¬ 
tions decrease with the order of the perturbation theory, 
as expected. 

We also study the scale uncertainties of both the FO 
and RESUM cross sections up to N 3 LO as a function of 
the center of mass energy y/s of the incoming protons 
at the LHC and our results are given in fig[3} Here, we 
vary hr in the range [O.lmn, lOmn] fixing hf = mu- In 
general, the scale uncertainties in both FO and RESUM 
results are found to increase with y/s precisely because 
of the increase in gluon fluxes. Irrespective of the order 
of the perturbation theory, the RESUM results are found 
to decrease the scale uncertainties remarkably compared 


to the FO results. Here, at N 3 LO, the cross sections will 
increase from hr = O.lmn to about hr = 0.5mn ( shown 
as solid lines in the Fig[3[ the dashed line corresponds to 
the one at hr = lOmn) and then start decreasing with 
further hr variation. Also for hr > mu, the N 3 LO cross 
section will decrease. Consequently for N 3 LO, the cross 
sections at the end points of the HR variation i.e. O.lmn 
and lOmn, will both be below the one at hr = wh. 



FIG. 3: Dependence of scale uncertainties in both the fixed 
order and resummed cross sections on y/s (see text). 


In conclusion, we have investigated the dependence of 
both the fixed order as well as the resummed predictions 
on the renormalization scale, using the recently available 
results on the Higgs boson production to N 3 LO in gluon 
fusion. For the resummed results, we systematically in¬ 
clude all the RG accessible logarithms, Lr, to all orders 
in the perturbation theory. While the fixed order N 3 LO 
result shows impressive scale reduction for the canonical 
choice of the renormalization scale between mn/2 and 2 
mu, there is still a significant dependence on the scale 
through these large logarithms which can spoil the be¬ 
havior if the renormalization scale is varied further away 
from this range. On the other hand, the resummed re¬ 
sults obtained in this letter show little dependence on 
the scale choice. For hr in the range [O.lmn, lOmn], the 
RG improved cross sections bring the scale uncertainties 
from about 27% down to about 1.5% at N 3 LO level. This 
approach can also be used for other processes such as top 
pair production, multi-jet production etc. 
Acknowledgments : We thank Nandadevi cluster com¬ 
puting facility at the Institute of Mathematical Sciences 
(IMSc) where the computation was carried out. GD 
thanks for the hospitality provided by IMSc where part 
of the work was done. GD also thanks P. Mathews for 
useful discussions and for his encouragement. 




























































5 


[1] C. Anastasiou, C. Duhr, F. Dulat, F. Herzog and 
B. Mistlberger, arXiv: 1503.06056 [hep-pli], 

[2] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716, 
1 (2012) ; S. Chatrchyan et al. [CMS Collaboration], 
Phys. Lett. B 716, 30 (2012). 

[3] H. M. Georgi, S. L. Glashow, M. E. Machacek, and 
D. V. Nanopoulos, Phys. Rev. Lett. 40 (1978) 692 ; 
A. Djouadi, M. Spira, and P. M. Zerwas, Phys. Lett. 
B 264 (1991) 440 ; S. Dawson, Nucl. Phys. B359 (1991) 
283 ; M. Spira, A. Djouadi, D. Graudenz, and P. M. Zer¬ 
was, Nucl. Phys. B453 (1995) 17 ; S. Catani, D. de Flo- 
rian and M. Grazzini, JHEP 0105 (2001) 025 ; R. V. Har- 
lander and W. B. Kilgore, Phys. Rev. D 64 (2001) 013015 
, Phys. Rev. Lett. 88 (2002) 201801 ; C. Anastasiou and 
K. Melnikov, Nucl. Phys. B646 (2002) 220 ; V. Ravin- 
dran, J. Smith, and W. L. van Neerven, Nucl. Phys. 
B665 (2003) 325. 

[4] S. Catani, D. de Florian, M. Grazzini, and P. Nason, 
JHEP 0307 (2003) 028. 

[5] U. Aglietti, R. Bonciani, G. Degrassi and A. Vicini, Phys. 
Lett. B 595 (2004) 432 ; S. Actis, G. Passarino, C. Sturm 
and S. Uccirati, Phys. Lett. B 670 (2008) 12 . 

[6] C. Anastasiou, C. Duhr, F. Dulat, E. Furlan, 
T. Gehrmann, F. Herzog and B. Mistlberger, Phys. Lett. 


B 737 (2014) 325 . 

[7] T. Ahmed, M. Mahakhud, N. Rana and V. Ravindran, 
Phys. Rev. Lett. 113 (2014) 11, 112002 ; T. Ahmed, 
M. K. Mandal, N. Rana and V. Ravindran, Phys. Rev. 
Lett. 113 (2014) 212003 ; S. Catani, L. Cieri, D. de Flo¬ 
rian, G. Ferrera and M. Grazzini, Nucl. Phys. B 888 
(2014) 75 ; T. Ahmed, N. Rana and V. Ravindran, JHEP 
1410 (2014) 139 ; T. Ahmed, M. K. Mandal, N. Rana and 
V. Ravindran, .JHEP 1502 (2015) 131 ; M. C. Kumar, 
M. K. Mandal and V. Ravindran, JHEP 1503 (2015) 
037. 

[8] J. Baglio and A. Djouadi, JHEP 1010 (2010) 064 . 

[9] V. Ahrens, T. Becher, M. Neubert and L. L. Yang, Phys. 
Rev. D 79 (2009) 033013 . 

[10] V. Ahrens, T. Becher, M. Neubert and L. L. Yang, Eur. 
Phys. J. C 62 (2009) 333 . 

[11] M. R. Ahmady, F. A. Chishtie, V. Elias, A. H. Fari- 
borz, N. Fattahi, D. G. C. McKeon, T. N. Sherry and 
T. G. Steele, Phys. Rev. D 66 (2002) 014010 . 

[12] S. Moch, J. A. M. Vermaseren and A. Vogt, Nucl. Phys. 
B 726 (2005) 317 . 

[13] A. D. Martin, W. .J. Stirling, R. S. Thorne and G. Watt, 
Eur. Phys. J. C 63 (2009) 189 . 

] 14] M. R. Whalley, D. Bourilkov and R. C. Group, hep- 
ph/0508110. 



